Intrinsic noise-induced phase transitions: beyond the noise interpretation. 
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We discuss intrinsic noise effects in stochastic multiplicative-noise partial differential equations, 
which are qualitatively independent of the noise interpretation (Ito vs. Stratonovich) , in particular 
in the context of noise-induced ordering phase transitions. We study a model which, contrary to 
all cases known so far, exhibits such ordering transitions when the noise is interpreted not only 
according to Stratonovich, but also to Ito. The main feature of this model is the absence of a linear 
instability at the transition point. The dynamical properties of the resulting noise-induced growth 
processes are studied and compared in the two interpretations and with a reference Ginzburg-Landau 
type model. A detailed discussion of new numerical algorithms used in both interpretations is also 
presented. 
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I. INTRODUCTION 

An important feature of nonlinear systems is their abil- 
ity to sustain an organized behavior even in the presence 
of a substantial amount of randomly fluctuating influ- 
ences. Even more strikingly, systems which in the ab- 
sence of fluctuations exhibit a disordered behavior can 
experience, under certain conditions, the emergence of 
spatiotemporal order upon addition of a suitable amount 
of noise [1]. The most basic manifestation of this fact 
is the existence of ordering phase transitions induced by 
noise in dynamical systems with spatial degrees of free- 
dom [2,3]. These transitions bring the system from a dis- 
ordered to an ordered phase as the intensity of the noise 
increases, contrary to naive intuition. By disordered (or- 
dered) phase we mean for example the homogeneous zero 
(non-zero) state corresponding to the coarse graining of 
a spin field with random (uniform) orientation. 

Ordering phase transitions are usually driven by multi- 
plicative noise terms, which depend on the system's vari- 
ables [4]. But the stochastic integrals associated with 
stochastic differential equations with multiplicative noise 
are not uniquely defined [5]. Among the many inter- 
pretations that can be given to these integrals, two are 
frequently used: the Stratonovich interpretation which 
follows the standard rules of calculus, but gives rise to 
non-intuitive statistical properties of the noise terms, and 
the Ito interpretation which avoids these problems, but 
at the expense of requiring new rules of calculus. 

Beyond the technical mathematical definitions, the 
physical implications of both noise prescriptions boil 
down to an important fact. The Stratonovich prescrip- 
tion for white noise yields the result one would get for 
a time-correlated noise in the limit of vanishing correla- 
tion time. The key point is that, as soon as the noise 
is slightly correlated, the stochastic variables defined by 
the corresponding Langevin equation build up correla- 
tions with the noise variable at equal time. This imme- 



diately implies that the multiplicative noise terms in the 
equation have a nonzero mean, even with a zero-mean 
noise. The result is the so-called Stratonovich drift, a 
net force induced by noise which is at the heart of most 
noise-induced phenomena, in particular concerning noise- 
induced ordering transitions. 

As for a given stochastic differential equation with mul- 
tiplicative noise the results do depend on the interpreta- 
tion, a preliminary analysis of the physical problem has 
to be performed to make a judicious choice. Our experi- 
ence indicates that there are a minimum of three possible 
situations: 

• If we start with a well established deterministic dif- 
ferential equation and some controlled parameter is 
allowed to fluctuate (experimental or realistic ex- 
ternal noise) one would always expect the noise to 
have a high-frequency cut-off and as a consequence 
the Stratonovich interpretation is usually argued to 
be the reasonable choice. 

• If the starting scheme is a master equation which 
is approximated by a Fokker-Planck equation, then 
one can write a stochastic differential equation with 
multiplicative noise in the Ito interpretation. This 
happens, for instance, in front propagation prob- 
lems on a lattice [6] . 

• Moreover quite often our initial scheme is a set of 
stochastic differential equations, and we would like 
to simplify the problem eliminating the most irrel- 
evant fast variables (those with a very short time 
scale). The interpretation of the flnal stochastic 
differential equation will depend on the order in 
which this procedure is performed with respect to 
the white noise limit. This is indeed a nontrivial 
task. 

Since the Stratonovich drift can drastically modify the 
behavior of systems, and since it may not always be 
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obvious what the appropriate noise prescription is in a 
given problem, it is particularly important to distinguish 
which noise effects are intrinsic, in the sense of occur- 
ring regardless the noise interpretation, and which ones 
are strictly associated to the Stratonovich drift. In other 
words, it is important to elucidate when the noise inter- 
pretation may only affect the quantitative behavior, and 
when it may indeed change the problem at a qualitative 
level. 

For the case of noise-induced phase-transitions, the 
noise prescription used so far in the literature is that 
of Stratonovich. Nevertheless, it could be argued that 
if the noise has an internal origin, one should in princi- 
ple expect Ito noise too, so it would be good to estab- 
lish whether, in the latter case, noise-induced transitions 
can occur. We will see that this is indeed the case for 
a recently discovered class of noise-induced phase tran- 
sitions. From a theoretical point of view, it is also im- 
portant to deal with Ito noise since then the continuum 
white-noise limit is either well defined or less singular 
than in the Stratonovich case [7]. This has important 
consequences in order to establish when the macroscopic 
observables will carry out a nontrivial, singular depen- 
dence on the spatial cut-off of the noise (Stratonovich 
case) and when such residual dependence will be weak 
(Ito case) [8]. 

Few contributions have appeared in the physics litera- 
ture on Ito calculus in extended systems. A comparative 
discussion about the mathematical problems involved in 
the two interpretations appeared in [9]. The role of the 
multiplicative noise in the Ito interpretation has been 
analyzed in the context of spatiotemporal intermittency 
[10] and front propagation [11]. Dynamical renormal- 
ization group calculations were presented in [12]. How- 
ever, noise-induced ordering phase transitions had been 
reported so far only in the framework of the Stratonovich 
interpretation [1-3]. In that case, the mechanism under- 
neath these transitions is that the multiplicative noise 
term has a non-zero average value, which produces a 
short-time instability of the disordered phase and induces 
the ordered phase to arise [2,3,13]. The instability can be 
linear [2,14] or nonlinear [3,15], but is in any case induced 
by the so-called Stratonovich shift. Due to the absence 
of such a drift, the Ito interpretation does not present 
this type of noise-induced ordered phase, or any other 
spatially ordered state [16]. 

Recently, however, a new type of noise-induced phase 
transition has been found which docs not occur via an 
instability of the disordered phase [17]. Here, the or- 
dered phase arises due to the balance between the re- 
laxing deterministic forces pushing the system toward 
the disordered state, and the activating multiplicative 
fluctuations pulling the field away from that state, in a 
type of entropy-driven phase transition (EDPT). This be- 
havior is the spatio-temporal extension of noise-induced 
transitions in purely temporal, zero-dimensional systems, 
where the probability distribution of the time-dependent 
variable exhibits a change in the number and type of 



its cxtrema as noise intensity varies [18]. A key idea 
in the model here studied is that the bimodality in the 
stationary probability density is not associated to a po- 
tential barrier, but has a dynamical origin. In fact, the 
dependence of the multiplicative noise term on the field 
is such that, for sufficiently large noise strength, the sys- 
tem escapes more easily from the central region than from 
the sides, despite the fact that the deterministic force al- 
ways drives the system towards the center. As a result 
the peaks of the probability density are off-center. An 
important difference with the usual bimodality associ- 
ated to a potential barrier is that in our case the char- 
acteristic relaxation time scales for the zero-dimensional 
model are of order one (0(£°)) as opposed to 0(exp(l/£)) 
which is characteristic of activation processes, s being 
a generic measure of the noise strength. In the spa- 
tially extended case, the spatial (diffusive) coupling of the 
field introduces an additional crucial ingredient, namely 
it freezes the domains impeding the fast relaxation pro- 
cess of the zcro-dimcnsional case. This gives rise to a 
well-defined, stable interface which then drives the much 
slower domain-growth dynamics. Since no Stratonovich 
drift is required to induce this effect (as opposed for in- 
stance to the case of Ref. [1] where it takes the form of an 
effective barrier) it is to be expected that the correspond- 
ing class of model exhibiting this behavior should also 
display noise-induced ordering in the Ito interpretation. 
In this paper we show that this is indeed the case, by 
comparing the behavior of the model introduced in [17] 
for both the Ito and Stratonovich interpretations with 
that of a standard Ginzburg-Landau model with mul- 
tiplicative noise (section II). We also analyze in detail 
the dynamical properties of the growth processes aris- 
ing from the noise-induced ordering transitions in the 
two cases (section HI) , which will be shown to share uni- 
versal characteristics (i.e. growth exponents) but differ 
in non-universal features (such as power-law prefactors) . 
Finally, algorithms that have been specially developed for 
generating the results presented in this paper, for both 
the Stratonovich and Ito interpretations, are described in 
detail in the Appendix. 

II. THEORETICAL ANALYSIS 

We will use a model of a class of systems for which 
the steady-state probability distribution can be obtained 
exactly. As a consequence, the existence of a phase tran- 
sition in this kind of systems can be studied without any 
dynamical reference. 

Our model corresponds to a relaxational flow in a 
free energy potential JF({0}), with a ficld-dcpendcnt ki- 
netic coefficient T((f)) and a fluctuating term fulfilling a 
fluctuation-dissipation relation [17]. The model is defined 
by the following stochastic partial differential equation: 
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We suppose that the noise ^{x, t) is Gaussian, with zero 
mean and correlation 



{^{x,t)^{x' ,t')) = 2a^S{x- x')6{t - t') 



(2) 



where cr^ is the noise intensity. Moreover, we choose the 
foUowing form for the free-energy potential T , 



j d^xi^Vo{^{S,t)) + ^ \Vcf>{x,t)' 



4d 



(3) 



Since we are dealing with spatially uncorrelated noise, 
we pcirform the analysis in a discrete space in order to 
avoid singularities [9]. In a d-dimensional square lattice 
of mesh size Ax and N = L''- cells, our model reads 



dt 'd4>i ^ 



(4) 



where only one index is used to label the cells, (f)i = (f){xi 
Ti = V[(j)i), and the noise satisfies the correlation: 



In discrete space, the free energy has the form 

D 



(5) 
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where the gradient term is approximated by the sum 

over nearest neighbors on the lattice in a standard way, 

IVf/'P Y.jenn+{i) ^'^Ax^^ , and nn+{i) stands for the 
d-nearest neighbors of i in the positive direction of each 
axis. For simplicity, we choose a monostable local poten- 
tial, 
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where a > 0. Finally, the kinetic coefficient r(^) is taken 
to depend on the field in the following way [17] 



m = 



1 + C4 



(8) 



This functional dependence of the kinetic coefficient fa- 
vors diffusion due to fluctuations in the disordered state. 

Our objective now is to study the equation (4) in the 
Stratonovich and ltd stochastic interpretations. The cor- 
responding Fokker-Planck equation for the probability 
density of the field t) can be written in a unified 

notation for both interpretations [5], 
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where i? = 1 for the Stratonovich interpretation and 
B = 2 in the ltd case. 

If no probability flux is present, the stationary solution 
Pst of (9) satisfies 

(dF Ba' dlnTA . _ . 

\dcl>i ^ 2Ax<i d(j>i ) ^* ^ Aar<i dcj^i ~ ^ ^ 

The solution of this equation is 

nt(W)~e-^^«^-''/-% (11) 

where we have introduced the effective free energy 
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Feff ({0}) = + ^ ^ In r, . (12) 

The above expressions can be written in continuum 
space as 



P,t({0})-e-^=«/-' 



(13) 



^•5^ff(W) ^ ^(W) + j d-^xlnVim) , (14) 



where (Tq^ = a"^ / Ax'^ stands for the effective noise inten- 
sity of a spatially white noise in a discrete space. 

We have thus seen that the stationary multivariate 
probability distribution can be obtained exactly in both 
the Ito and Stratonovich interpretations for the spatially 
extended EDPT model, and that both lead to very sim- 
ilar qualitative results. The only difference is an ex- 
tra factor 2 in the new term of the effective potential 
in the Ito interpretation. As is already known [17], 
the EDPT model presents a continuous ordering noise- 
induced phase transition in the Stratonovich interpreta- 
tion. But according to the results shown above, and as 
will be shown in the following Section, this model also 
exhibits an ordering transition in the ltd interpretation, 
although the location of the critical point will be differ- 
ent. We should remark here that, as in the case of the 
Stratonovich interpretation [17], this phase transition is 
not due to a short-time instability of the homogeneous 
null phase. Indeed, the linear equation for the first sta- 
tistical moment (^) can be computed to be [1] 



m 
dt 



[a + (2 ■ 



B)alc]{^) + ^^W^^). (15) 



For a > 0, the homogeneous null solution of this equa- 
tion is stable for all noise intensities, both for B = 1 and 
B = 2. Therefore, the mechanism of this phase transition 
must be different from the standard one. 
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III. STEADY-STATE BEHAVIOR 



A standard way of determining the existence of a noise- 
induced phase transition is by applying a mean-field ap- 
proximation to the Langevin or Fokker-Planck equations 
of the system [1,2]. In the present case, however, since 
we have obtained the exact multivariate probability dis- 
tribution in both interpretations, we will implement that 
approximation directly on the effective potential derived 
from (12). 

The mean-field approximation consists in replacing the 
exact value of the neighbor field in the Langevin or 
Fokker-Planck equation by a common mean-field value 
{<f>). In the present case, wc make such an identification 
in the neighboring values of the gradient term appearing 
in the effective free energy [see Eqs. (6) and (12)]: 
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In this way, the effective free energy becomes 



Feff({.^} , (</.)) = Vo{cj>,) + lnr(0O 

i=l ^ 



D 



N 



(.^i-(<^))n^^y,ff(^,,(^)). (17) 
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The unknown mean field value {(j)) is obtained self- 
consistently, according to 



/oo 
-oo 



(18) 



where the one-site probability distribution {Pst{{(j)}) = 



nili Psi{(l>i)) is given by 



-Veft/o-o^ 



(19) 



The mean-field predictions for {(j)) in the two-dimensional 
case are plotted in Fig. 1, where lines separating the sit- 
uations where (</)) = (disorder) and ((/)) ^ (order) 
are plotted for both the Ito and Stratonovich interpre- 
tations in the space of parameters D and ct^. The fig- 
ure shows that both interpretations predict a continu- 
ous noise-induced ordering phase transition, which oc- 
curs earlier (i.e., for lower noise intensities) in the Ito 
case. In particular, in the large coupling limit {D oo) 
the transition in the ltd interpretation takes place at a 
critical noise intensity (u^ = a/ Be, for Ao; = 1), that is 
half the critical value in the Stratonovich case, both of 
which coincide with the transition point in zero dimen- 
sional systems (with noise intensity a'^/Ax^) [18]. 
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FIG. 1. Phase diagram of the EDPT model, obtained 
from a mean-field analysis, in the Ito (continuous line) and 
Stratonovich (dashed line) interpretations. The horizontal 

dotted Unc corresponds to the value of D used in Fig. 2. The 
parameter values are a = 1, c = 0.5, and Ao; = 1. 

Note that, in contrast with the usual noise-induced 
transitions (which exhibit reentrant phenomena), the 
transition lines of Fig. 1 decay monotonously with 
(T^. This implies therefore that no minimum coupling 
strength is required in these models for a phase transi- 
tion to occur. 

In order to validate the results obtained from the 
mean-field approximation, we have performed extensive 
numerical simulations of model (4)- (8) in both the Ito 
and Stratonovich interpretations. To that end, we have 
developed a new type of numerical algorithm suitable for 
the implementation of both stochastic interpretations of 
the multiplicative noise. The derivation of this algorithm 
and a comparison with the well-known Heun algorithm 
(derived only for the Stratonovich interpretation) is pre- 
sented in the Appendix. The simulations have been per- 
formed on a square lattice of 256 x 256 cells of mesh size 
Ax = 1, with a time step At = 0.01 and periodic bound- 
ary conditions (except when explicitly indicated). Where 
necessary, wc have averaged over 10 realizations of the 
noise and the initial random conditions, corresponding 
to Gaussian or uniform distributions. In order to com- 
pute the mean field, we first evaluate the spatial average 
of the system: 



N 



(20) 



where N is the number of lattice cells, and (f>i{t) is the 
field value at the i cell. Once the spatial average reaches 
a stationary state, the temporal average is evaluated as 



Tm — Trr, 



Tm 



(21) 



t=T„ 



where Tm and Tm delimit the time interval within the 
steady-state regime in which the temporal average is cal- 
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culated. Afterwards, the realization average can be com- 
puted. 

The numerical simulation results for the two interpre- 
tations arc shown in Fig. 2, where they arc also com- 
pared with the predictions coming from the mean-field 
approximation. Due to the value of D chosen, the agree- 
ment between the mean-field estimate and the simula- 
tions is better for the Ito interpretation. In any case, the 
model exhibits a noise-induced ordering phase transition 
for both interpretations, as predicted by the mean-field 
approach. 
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FIG. 2. Mean-field and numerical simulation results for the 
EDPT model in the Ito (continuous line) and Stratonovich 
(dashed line) interpretations. Simulations have been per- 
formed for different system sizes: L = 16 (circles), i = 24 
(squares) and L = 32 (triangles) for Ito, and L = 64 (tri- 
angles), L = 48 (diamonds), L = 32 (squares) and L = 16 
(circles) for Stratonovich. D = A, and the rest of parameter 
values are those of the previous figure. 

We stress at this point that standard models exhibit- 
ing noise-induced phase transitions caused by short-term 
instabilities of the disordered phase do so only in the case 
of the Stratonovich interpretation. In order to illustrate 
this point, we present here for comparison what happens 
in the well-known case of the Ginzburg-Landau model 
with external multiplicative fluctuations [3]: 



= aq 



+ DV^(j) + (l)^{x,t) + r]{x,t) (22) 



where rj{x,t) and ^{x,t) are Gaussian and white noises. 
This system presents a noise-induced phase transition if 
we interpret the noise in the Stratonovich sense, but not 
if one uses the Ito interpretation. This can be seen in 
Fig. 3, where the two simulations share the same con- 
ditions and parameters. In the ltd interpretation the 
ordered parameter ((/>) remains always in the disordered 
state, due to the fact that the noise-dependent drift that 
causes the short-time instability is only present in the the 
Stratonovich prescription [16]. 
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FIG. 3. Bifurcation diagram of the Ginzburg-Landau 
model in the Ito and Stratonovich interpretations. Both 

mean-field and 2-d simulation results are shown. The pa- 
rameters used are L = 30, a = —0.2, & = 1, D = 4 and the 
additive noise intensity 0.5. 



IV. DOMAIN GROWTH DYNAMICS 



We have seen that the EDPT model in the presence of 
external fluctuations can reach a stationary ordered state 
described by a non-zero order parameter ((6) , for both the 
Ito and Stratonovich interpretations. This means that, if 
the system is initially in a disordered steady state {(f)) = 
corresponding to a small noise intensity, as the inten- 
sity of external fluctuations is increased above its critical 
value the system develops domains of the two new sym- 
metric stationary ordered phases, that grow with time as 
shown in Fig. 4. The figure shows that the system be- 
haves differently in the two stochastic interpretations for 
the same noise intensity, the Ito case being much more 
contrasted due to the fact that the order parameter is 
larger than in the Stratonovich case, which is very noisy. 

In this section we are concerned with the growth of 
these noise-induced domains. Although the mechanism 
that induces the phase transition is different from those 
that have been reported before, we can expect that, once 
the domains have appeared, their dynamics has the same 
characteristics as those of the domain growth following 
the quench of a system below its order-disorder transition 
temperature, as happens in the Ginzburg-Landau model 
[19]. 
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FIG. 4. Snapshots of evolving noise-induced domains for 
the EDPT model at t = 750 (left figures) and t = 1750 (right 
figures) in the Ito (top) and Stratonovich (bottom) interpre- 
tations. Parameters are: a = 1, c = 3, = 3.5 and L = 256. 

For non-conscrvcd order parameter models, one of the 
domains grows until it fills the whole system. The mecha- 
nism underlying domain growth in this case is the motion 
of the interface between domains caused by the inter- 
face structure. The translational velocity of the domain 
boundary has been found to be proportional to the mean 
curvature of the boundary, and independent of the free 
energy of the interface. This can be quantified by the 
equation of motion obeyed by the characteristic length 
(i.e. the average radius) of the domains of equilibrium 
phases, R{t), [20] 



_ r 



(23) 



where ^ is a model-dependent constant and F is the 
kinetic coefficient multiplying the diffusion term. This 
expression leads in a straightforward way to the AUen- 
Cahn law of domain growth: 



R{t) oc V2Ar t^'"^ 



(24) 



In the time regime where this law is verified, R{t) is the 
only characteristic length of the system, and a scaling be- 
havior for its spatial structure at different times is found. 
All these results are known to apply also in the case of 
standard noise-induced phase transitions caused by linear 
instabilities of the homogeneous disordered phase [14]. 
Wc want to find out whether the same thing happens in 
the EDPT described in this paper. 

In order to characterize the dynamics of model (4), 
wc let our system evolve from an initial disordered state, 
and compute the isotropic correlation G(r, t) function at 
different times. We use the following normalization 



G{r,t) 
G(0,i)- 



(25) 



Let us consider a time regime in which there is only one 
characteristic length R{t) in the system, which is related 
to the average size of the domains. There arc several 
possible definitions for i?(t), but all of them should lead 
to the same results. We have chosen R{t) as the dis- 
tance at which g{r, t) has half its maximum value. In 
this time regime, we can apply the scaling hypothesis for 
a d-dimensional system 



g{r,t)=g{r/R{t)), 



(26) 



with no other explicit time dependence. When these re- 
lations hold, the spatial structure of the system at dif- 
ferent times is statistically equivalent, except for a scale 
factor. Since the domain growth is more clearly observed 
far from the critical point, we have taken new param- 
eter values accordingly. The numerical results in the 
Stratonovich interpretation for the scaled pair correlation 
function are represented in Fig. 5. As it was shown in 
[19], the pair correlation function exhibits a discontinuity 
in its first derivative in the presence of noise sources. We 
have eliminated this discontinuity by fitting a parabolic 
function in the origin (r = 0, Ax). The same study has 
been made in the case of Ito interpretation under the 
same conditions and parameters. 
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model in 



Scaled pair correlation function for the EDPT 
the Stratonovich interpretation for t=1300 (cir- 
cles) and t=2000 (squares), and in the Ito interpretation for 
t = 1200 (triangles) and t = 1800 (diamonds). The parameter 



values are o = 1, c = 3, Z) = 4, 



3.5 and Aa; ; 



We now compare the temporal evolution of the char- 
acteristic length of the system R{t) for the two stochas- 
tic interpretations. Fig. 6 presents this comparison for 
equal values of the noise intensity. From these numeri- 
cal results we can conclude that the Allen-Cahn law is 
satisfied for the two interpretations, and that there is 
a time regime in which the system is self-similar. One 
interesting fact is that domain evolution in Ito is slower 
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than in Stratonovich. and in both cases much slower than 
the Ginzburg- Landau model. This fact can be explained 
looking at the constant prefactor ^2 AT of the AUen- 
Cahn law (24). In the Ginzburg-Landau model F = 1, 
but in the EDPT model this quantity is field dependent 
(8), and can be approximated by 



(27) 



1 + C(<^2^ 1 + C(<^)2 ■ 



According to this expression, and since for a fixed (7^ we 
have that (</>)/ > (0)s, as a consequence we should ex- 
pect the slowest growth for the Ito EDPT case, and the 
fastest one for the Ginzburg-Landau model. This is what 
we can see in Fig. 6. 
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FIG. 6. Allen-Cahn law for the GL model (same parameter 
values as in Fig.3 with e = 10"'' and cr^ = 0.6) and the EDPT 
model for equal noise intensities (and different mean fields). 
The latter is computed in both the Ito and Stratonovich in- 
terpretations. The parameter values for the EDPT model are 
o = 1, c = 3, D = 4, = 3.5 and Aa; = 1. 



In order to eliminate the influence of the stationary 
mean-field value {(f)) on the growth rate, we have com- 
pared the evolution of the system under the two inter- 
pretations using in each case a different noise intensity, 
so that the mean field has the same value in the two 
cases. The results are shown in Fig 7, where we have 
fixed {(j)) = 3.15 for the two interpretations, for which 
we need = 16 in the Stratonovich interpretation and 
cr^ = 6 in the Ito interpretation. As can be seen, in both 
interpretations the system seems to evolve at the same 
rhythm, although the slope in Ito is slightly higher than 
in Stratonovich. We think that this small difference is 
due to the fact that although {(j)) is the same for both 
interpretations, the Stratonovich case is much more fluc- 
tuating, and hence we have to expect a larger (t/)^) and 
accordingly a lower slope. 
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FIG. 7. Allen-Cahn law for the EDPT model for equal 
mean fields (and different noise intensities), under both the Ito 
and Stratonovich interpretations. For Stratonovich = 16, 
for Ito cr^ = 6. Other parameter values are: a = 1, c = 0.5. 



V. COMMENTS AND CONCLUSIONS 

It is worth commenting here that an effective model 
can be developed which has the same stationary solu- 
tions as the EDPT model, but different dynamics. The 
dynamical equation for this effective model with Aa; = 1 
is 



dt 



where, as before, B is a parameter whose value indicates 
the interpretation that we are mimicking (B = 1 for 
Stratonovich and B = 2 for Ito) . The correlation of noise 
is given by Eq. (2). The equation of motion of the mean 
value of the field in the linear approximation is 



^ = (Sa2c-a)(0) + ^V2(0), 



(29) 



which tells us that, for tr^ > a/ Be, the homogeneous 
phase {(j)) = is unstable. This instability does not ap- 
pear in the EDPT model (see Eq. 15). According to 
this result, we have to expect an initial transient faster 
in this model (as in the Ginzburg-Landau model) than in 
the EDPT cases. This fact has been checked numerically 
and it can be seen in Fig. 8. We can clearly appreciate 
that the effective model has a much faster initial tran- 
sient that the EDPT model, for the same values of the 
parameters which is a signature of the different character 
of the instability of the initial state. While this observa- 
tion applies also to the zero-dimensional version of the 
model, the crucial ingredient in our EDPT model is the 
role of the spatial coupling, which prevents the fast tran- 
sition between the two probability peaks. The fact that 
this transition occurs in a deterministic time scale in the 
zero-dimcnsional case is what distinguishes the problem 
from the usual, barrier-crossing bistability. In the spa- 
tially extended case, however, it is precisely the spatial 
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coupling what generates an effective barrier allowing for 
the formation of stable domains, with an interface-driven 
dynamics. 
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FIG. 

mn =< 



8. 



Transient evolution of the quantity 
■^{t) > / < 4> which measures the emergence 
of order from homogeneous initial condition c6 = 0. The let- 
ter I means the Ito case and the value of the intensity of the 
noise is inside the parenthesis. 

In conclusion, we have presented a nonequilibrium field 
model for which one can compute exactly the stationary 
probability distribution, and which exhibits an intrinsic 
noise-induced ordering phase transition irrespective of 
the stochastic interpretation of the multiplicative noise 
term. In particular, the phase transition is found in the 
Ito interpretation, where so far, noise had only been seen 
to have a disordering effect. The same model can be stud- 
ied changing the diffusive term by the spatial coupling of 
the Swift-Hohenberg model in which case a noise-induced 
pattern transition is found [21]. These type of models 
do constitute a generalization of the Horsthemke-Lefever 
noise induced transitions to genuine noise induced phase 
transitions in extended systems. 
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APPENDIX A: STOCHASTIC ALGORITHMS 

Here we will derive an alternative algorithm which is 
an extension of the well-known Heun algorithm, valid for 
both the Ito and Stratonovich interpretations of stochas- 
tic differential equations with multiplicative noise. Our 
aim is to simulate numerically the following stochastic 
differential equation on a d-dimensional lattice. 



d(j){x, t) 
dt 



f{cj>{x, t), V) + g{ct>{x, tmS, t). (Al) 



First of all, we write this equation in a discrete space as 
follows. 



d<i>i{t) 
dt 



(A2) 



where i stands for the position inside the lattice, and the 
noise correlation is given by Eq. (5). 

The first step in the derivation of the algorithm is to 
integrate formally Eq. (A2) to get 



/t+At 
hm'))dt' 

g{mW)dt' (A3) 



The first integral in (A3) is evaluated according to a 
second-order predictor-corrector algorithm. 



L ■ 



g{m)mdt' (A4) 



where 4>{t) is the predictor term defined as the first-order 
solution of (A3), 

Ut) = Mt) + fi{<i>i{t))^t + gim))Xi (A5) 

This expression defines the first equation of the algo- 
rithm, which is independent of the stochastic interpre- 
tation. Xi(t) is the Wiener process, defined as 



/t+M 
Ut')dt' 

and whose numerical implementation is 



/2a2At 



(A6) 



(A7) 



where on are independent Gaussian random numbers of 
zero mean and unity variance, and they are implemented 
using [22]. 

The second integral in Eq. (A3) is not well defined, 
and one needs to make a prescription for its evaluation, 
at least up to first order in At. 

The standard Heun algorithm works for the 
Stratonovich interpretation, and makes the following 
assumption 



(AS) 
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Accordingly, the second equation of this algorithm is 

(A9) 



On the other hand, in the Stratonovich calculus this in- 
tegral is interpreted as [5] 



/•t+At 

gm'mt')dt' = g 



(AlO) 



so that the second equation of the algorithm is 
*(, + At)=«t, + »))+^(Mi, 

^Jm^]m>. (All) 



which is not exactly the standard Hcun algorithm. This 
is the algorithm that has been used in this paper for 
the Stratonovich interpretation results. Its advantage is 
that, in contrast to Heun algorithm, our method has an 
analogue in the Ito interpretation, for which the same 
integral is defined as [5] 

/t+At 
9im')Mt')dt' = gi{cP{t))Xi{t). (A12) 

Therefore, the second equation of the algorithm in Ito 
interpretation reads 

+9,(0(t))X,(t). (A13) 

Given these results, the algorithm proceeds by evaluating 
first the predictor contribution (A5) and, using this value, 
computing the corrector term (A9), (All) or (A13), cor- 
responding to the Hcun, Stratonovich or Ito algorithms 
respectively. All these three different algorithms are ap- 
proximations up to the same order (second order in the 
deterministic part but first order in the stochastic one), 
when properly expanded in powers of At. One can check 
that there are no differences, up to these orders, between 
the Heun and Stratonovich algorithms, as it should be. 
Nevertheless, the Stratonovich prescription has an ex- 
tra term, l/2g((j)i)g' {(pi)Xi(t)'^ , with respect the Ito one, 
which is of order At. Our Ito algorithm also agrees with 
the one presented in [23] up to order At^. 
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